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Abstract 



Recently, Neubert has suggested that a certain class of nonperturbative cor- 
rections dominates the shape of the electron spectrum in the endpoint region 
of semileptonic B decay. Perturbative QCD corrections are important in 
the endpoint region. We study the effects of these corrections on Neubert's 
proposal. The connection between the endpoint of the electron spectrum in 
semileptonic B decay and the photon spectrum in b — > 57 is outlined. 
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I. INTRODUCTION 



The electron energy spectrum near its endpoint in semileptonic B meson decay arises 
from b — > u transitions and provides one method for the extraction of the Kobayashi- 
Maskawa mixing angle V u b from experiment. The spectrum must be known accurately within 
a few hundred MeV of its endpoint, since it is only in this region that the large background 
due to the dominant b — > c weak transition is kinematically forbidden. Thus, the separation 
of the rare b —>■ u decay from the inclusive spectrum relies upon a theoretical understanding 
of the shape of the spectrum in this small region. Unfortunately, it is precisely this region 
which is the least well understood theoretically. 

The endpoint region of inclusive semileptonic B decay has been studied extensively. The 
first approaches relied on QCD models. Grinstein et al. [jl] used a constituent quark model 
to sum over exclusive charmless final states in this region, assuming that the spectrum is 
dominated by a few low-lying resonances. Altarelli et al. |J computed the spectrum in the 
free b quark decay model, augmented by the inclusion of a model of the Fermi motion of the 
b quark in the B meson. More recently, a QCD-based approach has been formulated in the 
context of the heavy quark effective theory (HQET). Using an operator product expansion 
(OPE) and the HQET, Chay, Georgi and Grinstein || have shown that the free b quark 
decay model describes inclusive semileptonic B decay to leading and first subleading order 
in a systematic expansion in 1 /nib, where nib is the b quark mass of the HQET. The first 
non- vanishing corrections to the free quark decay result are of order 1 /ml, and have now 
been computed These corrections arise from higher order terms in the OPE whose 

matrix elements contain information about the state of the b quark inside the hadron. 

At leading order, the electron spectrum is governed by quark kinematics with an endpoint 
at E e = nib/ 2, rather than at the physical endpoint Mb/2 which is determined by the 
B meson mass Mb- The higher order terms in the 1 /nib expansion produce corrections 
to the free quark decay spectrum, causing it to "leak" beyond the free quark endpoint. 
Understanding this process is crucial for extracting V u b, since the difference (M B — nib)/ 2 is 
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expected to be several hundred MeV, and is comparable to the 330 MeV energy difference 
between the b — > u and b — > c endpoints. Recently, Neubert has shown that the most 
singular terms in the l/m& expansion can be used to define a "shape function" of the 
spectrum, which is determined by a certain set of nonperturbative matrix elements and is 
model-independent ]7]]. This shape function describes the electron energy spectrum beyond 
the kinematic endpoint of the free quark decay (neglecting QCD radiative corrections). 
In this paper we examine the influence of perturbative QCD corrections on the endpoint 
region. These corrections are particularly important due to the presence of a Sudakov 
double-logarithmic suppression of the free quark decay rate at the endpoint. 

This paper is organized as follows. In Section 2, we review the operator product expansion 
analysis of the differential decay width for the endpoint region, neglecting perturbative QCD 
radiative corrections. The summation of the leading nonperturbative singularities to the 
shape of the endpoint spectrum is presented. We show that this summation can be obtained 
from the free b quark decay result by suitably averaging the free quark decay result over the 
residual momentum of the b quark inside the B meson. Since the leading nonperturbative 
corrections can be generated by this procedure, radiative corrections can be included by 
computing radiative corrections to the free quark decay result and then averaging over the 
residual momentum of the b quark. In Section 3, we consider the radiative corrections to free 
quark decay and show how they modify the shape of the endpoint of the electron spectrum. 
Numerical results and conclusions are presented in Section 4. 

II. LEADING NONPERTURBATIVE SINGULARITIES 

The inclusive differential decay distribution for B — > X UjC e V is determined by the imag- 
inary part of the time-ordered product of two weak currents, 

= ~ij d 4 xe- iq - x (B\ T{J"\x), J w {0)} \B) , (2.1) 

where J M = g7 M (l — 7 5 )6 and q = u,c. The time-ordered product may be expanded in 
inverse powers of the b quark mass using an operator product expansion ||, and in powers 
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of a s (m&). In this section we will concentrate on the l/m b expansion. From the operator 
product expansion of the hadronic tensor, one obtains an expression for the inclusive electron 
energy spectrum, dT/dy, where y is the rescaled electron energy, y = 2E e /rrib. The leading 
term in the 1/rrib expansion produces the result of the free quark decay model, in which 
the inclusive semileptonic decay rate is given by the decay of a free, on-shell b quark. The 
endpoint of the electron spectrum is at y = 1. The sub leading terms represent corrections to 
free quark decay, in which certain features of the motion of the b quark inside the B meson 
are taken into account. The expansion is in powers of 

e = A/m b , (2.2) 

where A is a scale typical of the strong interactions of QCD, perhaps 300 to 500 MeV. 

Neglecting perturbative a s (m b ) corrections, the electron energy spectrum for B — > X u e V 
decay is given by 



If = {2(3 - 2yW + 4(3 - ytf* - ^±M Kb - Ml|±M Gt } , (I _ y) 

r 4 16 1 2 

+ \2E b - -K h + -G b j 6(1 —y) + -K b S'(l - y) , (2.3) 

up to corrections of order e 3 , where T is the free quark decay width 

and 9(x) is 1 if x > and zero otherwise.^ E b , K b and G b are hadronic matrix elements of 
order e 2 , defined by 



E b = G b + K b . 



D 
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K b = (B(v)\b v —^b v \B(v)) , 
Zm b 

G b = (B(v)\b v g^f^b v \B(v)) , (2.5) 
4m£ 



*Eq. ( |2.3| ) holds for massless leptons. Lepton mass effects may be included ||, but they do not 
change the behavior of the endpoint spectrum in an important way. 
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where b v is the b quark field in the HQET. The factor of 9(1 —y) in the first term is required 
because the tree level decay distribution does not vanish at the boundary of the Dalitz plot. 
The 5(1 — y) and 5'(1 — y) singularities arise because some higher order terms in the 1/rrib 
expansion have the form of derivatives with respect to y of lower order terms. Since the 
free quark decay distribution does not vanish at the endpoint, this generates singular terms 
in the decay spectrum. These singularities imply that the l/m& expansion breaks down at 
y=l. 

Eq. (|2.3| ) is the decay spectrum including all corrections of order l/m|. To all orders 
in l/m&, the decay spectrum dT/dy obtained from the OPE at zeroth order in a s has the 
structure 

ir^ = *(l-v) ( e ° + 0e + e 2 + ---)+5(l-y) (0e + e 2 + -.-)+5'(l-y) (e 2 + e 3 + ■ ■ •) 
+ • ■ ■ + 5^(1 - y) (e n+1 + e n+2 + ■■■)+•■•, (2.6) 

where e n denotes a term of that order, which may include a smooth function of y. It is a 
nontrivial prediction of the heavy quark effective theory that the terms proportional to e 
in this expansion vanish [|J, as is evident in eq. (|2.3|) . Although the theoretical expression 
for dT/dy is singular at the endpoint y = 1, the total semileptonic width is not. The 
contribution to the total rate of a term e m 5^ n \l — y) is of order e m , so the semileptonic 
width has a well-behaved expansion in powers of l/m&, 

r = r (l + Oe + e 2 + e 3 + . . .) , (2.7) 

where the term proportional to e vanishes. 

The semileptonic decay width for b — > u is difficult to measure because of background 
contamination from the dominant b —>■ c semileptonic decays. It is therefore important to be 
able to compute the semileptonic decay rate for b —>■ u transitions near the endpoint y — 1, 
since the kinematic endpoint of the b —>■ c spectrum is below the b — > u endpoint. One 
way to calculate the endpoint spectrum is to weight the differential distribution dT/dy by a 
normalized function of width a around y — 1. We will refer to this procedure as "smearing." 
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Most of the details of the smearing procedure are unimportant; the only quantity of relevance 
is the width a of the smearing region. A physically meaningful result can be obtained by 
smearing over a large enough region in y such that the singular corrections to dT/dy are 
small. In ref . Q , it was shown that the singular corrections are small if the smearing width 
is chosen so that o ^> e. We will now show that by summing the leading singularities, one 
can choose a of order e. 

The singular distribution e m 8^ n \l — y) (where m > n) smeared over a region of width 
a gives a contribution of order e m /a n+1 to dT/dy. If the width a of the smearing region is 
of order e p , the generic term e m 5 (jl \l — y) yields a contribution of order e m -( n + 1 )p. Since 
m > n, this shows that the l/m& expansion for the spectrum breaks down unless p < 1, i.e. 
the smearing region cannot be made narrower than of order e. If p > 1, the l/rnj, expansion 
breaks down because it is dominated by an infinite number of terms at large values of n. This 
divergence is not associated with the failure of the OPE due to the presence of resonances 
with masses of order the QCD scale ||. The region in which such resonances dominate the 
final state is of width e 2 , while the expansion breaks down upon smearing over any region 
of size e 1+s , where 5 > 0. 

If the smearing region is chosen to be of order e, the form of the expansion ( |2.6| ) shows 
that the leading terms of the form 9(1 — y) and e n+l 5 (jl \l — y) all contribute at order unity 
to dT/dy, all terms of the form e n+2 5^ n \l — y) contribute at order e, etc. Thus one can, in 
principle, obtain the decay spectrum smeared over a width of order e if one can sum the 
leading singularities in eq. ( |2.6|) . The sum of the leading singularities produces a distribu- 
tion dT/dy of width e, and with a height of the same magnitude as the free quark decay 
distribution for dT/dy, i.e. with a height of order one. The subleading singularities produce 
a distribution which is also of width e, but has a height of order e times the distribution 
obtained by summing the leading singularities. The decay distribution dT/dy cannot be 
obtained with a resolution finer than e without summing all the subleading singularities. 

Neubert has shown that the series of leading singularities 
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~ = A 6(l-y) + 0e 5(1 - y) + A 2 e 2 5'(1 - y) + • ■ ■ (2.8) 
i o ay 

may be resummed into a "shape function" , which describes the behavior of the theoretical 
spectrum in the region beyond the free quark decay endpoint at y = 1 0. These terms arise 
in a particularly simple way in the OPE, because they come only from the expansion of the 
quark propagator which connects the two currents. The shape function has a width of order 
e and height of order one. 

The series of leading singularities ( |2.8| ) can be obtained by averaging the free quark decay 
result over the residual momentum of the b quark in the B meson . This simple procedure 
is important since it will also enable us to obtain the leading nonperturbative singularities 
for the radiative corrections by only calculating radiative corrections to free quark decay. 

The differential decay distribution is obtained from the tensor T^ v defined in eq. Q2.1| ). 
This tensor is a function of the momentum transfer to the leptons, q, and the velocity of 
the B meson, v. The differential decay distribution is proportional to the hadronic tensor 
contracted with the lepton tensor L^, which depends on the electron and neutrino momenta, 
k e and k v : 

, « W^L, U , (2.9) 

ax ay dq A 

where W^ v is the discontinuity of T^ v across the physical cut, W^ v = —ImT^/ir. The 
constant of proportionality in eq. fl2.9p involves G 2 F and the mixing angle |Kifc| 2 . The dimen- 
sionless variables x, y and q 2 are defined by 

2k v -v 2k e -v ^ q 2 

x = , y= , q = —. (2.10) 

rrib mi, m t 

The lowest order (in l/m b ) decay distribution dT iree /dx dy dq 2 is the decay distribution for 
a free on-shell b quark with mass m& and the same velocity v as the B meson. However, 
the b quark in the B meson is off-shell with a distribution of residual momentum k. The 
off-shell b quark, with momentum m^v + k, may be viewed as an on-shell quark with mass 
m' b and velocity v', where m' b v' = m^v + k. The decay rate for such a quark is obtained by 



evaluating the lowest order expression for dT{ ree /dy in the rest frame of the moving quark, 
and then boosting back to the rest frame of the B meson, 

dT = -L rfr free (x', y', q'\ m' b ) . (2.11) 

Note that all scaled quantities depend implicitly on m&, and hence must be primed. We now 
replace m' b v' — > m^v + k and average over the residual momentum k^. Expanding in k^/nib, 
we obtain a series of the form Q 

/ [l + 2t»-A;/m;, + A; 2 /mg] 1/2 " 



(dr> 



1 + v ■ k/mb 



d 1 9 9 
1 + ^- + -^F 2 - + 

orribV^ 1 2 orribV^ 1 orribV^ 2 



dTf ree J 
(2.12) 



where (■) denotes an average with respect to the distribution of the momentum k of the b 
quark in the B meson. The derivatives with respect to m^v 11 can be rewritten as derivatives 
with respect to x, y and q 2 using the chain rule. Terms with n derivatives with respect to 
rribV^ in eq. (|2.12| ) turn into terms with n x , n y and n q derivatives with respect to x, y and q 2 
respectively, where n x + n y + n q < n. The expansion of dT/dy is then obtained by integrating 
the expansion of dT/dxdydq 2 with respect to x and q 2 . The explicit computations to order 
1 / m 2 are given in ref . Q . 

In this paper, we are interested in summing the most singular terms in dT/dy near y — 1 
to all orders in l/m&. These terms are found by retaining the terms in eq. ( |2.12| ) with 
the maximum number of y-derivatives at each order in l/m&. This corresponds to only 
retaining the d n /dy n term in d n / drribV^ 1 . . . drribV^" in eq. fl2.12| ) and ignoring the prefactor 



1 /2 

[1 + 2v ■ k/mb + k 2 /ml] / [1 + v ■ k/rrib]- Terms with derivatives with respect to x or q 2 do 
not generate derivatives with respect to y on integration over x and q 2 , and are less singular 
than the terms we have retained. The most singular terms are thus obtained using 



V ( dy \ " ( 2 , ; \ " „. : , / 2 , ; 



n 



(K^-yv^) (k e -v)^\ , (2.13) 



ydmbV^ J \dmbV^ dy ) \rrib dy J \rrib dy J 

which gives the leading singularities, 
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dT dT ivec , / 2 \ j. 9 /"dTfre 



+ 



% % VmJ e 1X1 dy \ dy J 

+ " ' + ' ' ' (£)" {k - •>« -(*•- fif) + " ' 

where fc e = k e /rrib. Eq. ( |2.14 ) sums the leading nonperturbative corrections in the endpoint 



region, provided one interprets the residual momentum k in eq. ( |2.9j ) as the operator iD and 
the average as the expectation value of the resulting operator in the 5-meson state. There 
is no operator ordering ambiguity for the leading singularity in this identification, because 
D^i . . . j s contracted with the completely symmetric tensor (k e — t>) m • • • (k e — v)^ n , 
and so the commutator D u ] does not contribute. Only the part of the matrix element 
(B(v)\iD^ ■ ■ ■ iD^ n \B(v)) proportional to the tensor structure t> Ml • • contributes to the 
most singular terms, since (k e — v ) 2 vanishes at y = 1 [^]. Neglecting perturbative a s (ra&) 
radiative corrections, the most singular terms in eq. ( [2.1 4|) are ^-functions and their deriva- 
tives, which arise from differentiating the factor of 6(1 — y) in dT^ ee /dy. Dropping the n = 
term in eq. (|2.14 ) and allowing the derivatives to act only on the ^-function gives Neubert's 
shape function 

oo ™ Qn 

S (v) = E ZHTTT A - «)w • • • & - (B{v)\iL»* ■ ■ ■ iD^ \B(v)) — 6(1 - y) . (2.15) 

n=l m b n - °y 

This procedure for averaging over residual momentum produces the same result for the 
leading singularities as the operator product expansion. As discussed in ref. [|J, one can 
use reparameterization invariance || to show that averaging over residual momentum gives 
the same answer as the OPE, provided one neglects the commutator [D**, D u \ and higher 
dimension operators involving light quark fields. The commutator and higher dimension 
operators do not contribute to the most singular terms, and so averaging over residual 
momentum will be adequate for this discussion. 

It is simple to understand how this averaging procedure generates a shape function which 
extends beyond the free quark decay endpoint. If the energy of the b quark is allowed to 



9 



fluctuate from its on-shell value, occasionally it will have an energy larger than its free value 
m&. This fluctuation corresponds to a situation in which the quark has temporarily absorbed 
some energy from the light degrees of freedom in the B meson; if it decays weakly at this 
moment, then an energy E e > rrib/2 may be given to the electron. 



III. RADIATIVE CORRECTIONS 

The advantage of the averaging procedure for obtaining the leading nonperturbative 
singularities as y — > 1 is that it generalizes straightforwardly to the case when radiative 
corrections are included. The averaging procedure applied to the free quark decay distribu- 
tion including radiative corrections yields the leading nonperturbative singularities including 
radiative corrections. 

The one-loop QCD contribution to the free quark decay process, including both virtual 
gluons and real gluon emission, has been computed |1(J. The corrected electron spectrum 
takes the form 

d^ivee dT 



dy dy 



' 2as -G(y,m q ) + 0(a 2 s ) 



(3.1) 



3tt 

where T is the tree level free quark decay rate. Perturbative QCD corrections do not extend 
the electron spectrum beyond the free quark decay endpoint y = 1. This can only occur 
because of the nonperturbative l/m^ corrections discussed in the preceding section. In the 
interesting case m q = relevant to the transition B — > X u eV, G(y, 0) is given by fT0] . |lT|f 



31 5 

G(y, 0) = G(y) = ln 2 (l - y) + — ln(l - y) + tt 2 + - + (vanishing as y -+ 1) . (3.2) 

6 4 

The leading singularity at each order in perturbation theory is proportional to a"ln 2n (l — 
y). These singularities lead to a breakdown of the perturbative QCD expansion near the 
endpoint y = 1, unless they can be summed. The double logarithms have been shown to 
exponentiate 03], yielding an expression which formally has the structure 



^The expression for G(y) is taken from ref. [11| 
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-^ = R(y)-°, 3.3 
dy dy 

where 

E( y ) = exp{-^ln 2 (l-y)} . (3.4) 

This is the form for the decay spectrum used by Altarelli et al. 0. The Sudakov form factor 
R(y) causes the electron spectrum to vanish at the free quark endpoint y = 1. 

The contribution to the endpoint shape of the electron energy spectrum coming from the 
exponentiated double-logarithm in R(y) is a calculable effect. One might hope that once 
this leading radiative correction has been accounted for, it would be consistent to include 
the leading higher dimension operators using eq. (|2.14 ) and neglect all sub leading radiative 



corrections. However, we find that for very large mj this is not the case; the perturbative 
expansion is so poorly behaved at large orders in a s (rrib) that it is necessary to sum an 
infinite number of infinite series before including nonperturbative effects with eq. ( ^.14 ). 



Nevertheless, for the case of interest m& w 4.5 GeV, neglecting the subleading radiative 
corrections may provide a reasonable approximation for the endpoint of the electron energy 
spectrum. 

Before analyzing the general structure of the radiative corrections, it is instructive to 
consider a simple example which illustrates the importance of subleading radiative correc- 
tions. Consider the order a s correction given in eq. ( |3.1|) . This correction has ln 2 (l — y) 
and m(l — y) singularities as y — > 1. The ln 2 (l — y) singularity is summed into the Sudakov 
form factor R(y), leaving the subleading ln(l — y) singularity. This subleading logarithmic 
singularity must also be understood in order to determine the effect of radiative corrections 
on the endpoint energy spectrum [13]. To see this, note that it is possible to write two differ- 



ent expressions for the decay spectrum which contain the same Sudakov leading singularity, 
but which have very different behaviors as y — > 1. The first expression is the conventional 
definition |2[] 



dTfr ee _ R(y\ 

dy dy 



3tt 
11 



G{y) 



(3.5) 



where 



(3.6) 



However, one can also rewrite the decay spectrum as 



dy 



dy . 



In ~ 

R{y) - ^G{v) 



(3.7) 



which is equally valid to order a s . The two expressions (pTq) and (pTq ) have the same 
ln 2 (l — y) singularity as y — > 1, but differ in the subleading terms. The first expression 
(p.5|) vanishes as y — > 1, whereas eq. (|3.6| ) diverges as ?/ — ► 1. Thus, the exact form of the 
subleading singularity is required in order to determine the shape of the spectrum very near 
the endpoint. 

We will now demonstrate that in the limit m& — > oo, summing the most singular l/m& 
corrections with eq. (|2.14 ) cannot be used to improve the behavior of the electron spectrum 
near y = 1 without first summing an infinite number of subleading perturbative QCD 
singularities. In a schematic notation in which we include only the powers of a s and ln(l— y), 
the radiative corrections near y = 1 have the structure 



+ a s ln 2 (l - y) + a s ln(l - y) + a s 

+ a\ ln 4 (l -y) + a 2 s ln 3 (l - y) + a 2 s ln 2 (l - y) + a\ ln(l - y) + a 2 s 
+ a 3 s ln 6 (l -y) + a 3 s ln 5 (l - y) + a 3 s ln 4 (l - y) + a 3 s ln 3 (l - y) + ■ ■ 



+ 



(3.8) 



The first column, containing terms of the form a" ln 2n (l — y), exponentiates into the Sudakov 
factor R(y), after which the most singular terms remaining are of order a n s In 2 "' 1 (1 -y). We 
may write the mth column of the expansion ( |3.8| ) as an infinite series of the form 

(3.9) 



C m (y)= E WCln 2 "-™ +1 {i-y). 



n=[m/2] 
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The series of leading singularities corresponds to m = 1; for this case, and only this case, 
the coefficients 

1 / 2 \ n 

h m = - ] — (3.10) 



n! V 37i\ 

have been computed for all n, and the sum Ci(y) is The series C m (y) for m > 1 

represent an infinite set of infinite series, for which the behavior of the coefficients b mn for 
large n is not known. 

The unknown sub leading series C m (y), m > 1, limit the accuracy with which one can 
determine the electron energy spectrum. For perturbation theory to be valid, one has to 
remain in a region in which all the subleading terms are small, since their structure is not 
known, i.e. all the terms beyond the first column of eq. (|3.8|) must be small. This condition 
requires that a"ln 2n_m (l — y) <C 1 for all n and all m > 1, or that a s C 1 and 



a s \n\l-y) < 1, (3.11) 



which is the condition required for n — > oo with m fixed. If eq. ( |3 . 1 1| ) is satisfied, the 



first column sums to R(y), the second column is of order ^foT s times the first column, the 
third column is of order ^faT s times the second column, and so on. The condition (3.11) 



has converted the QCD perturbation series in eq. ( |3.8| ) into an expansion in yfoTg. Summing 
all the leading singularities a"ln 2ra (l — y), or summing any finite number of columns of 
eq. (|3.8|) , does not increase the region of validity of the perturbation expansion, since the 
condition that the next column be small is still eq. ( |3.11|) . To increase the region of validity 
of the perturbative expansion, one must sum all the terms of the form a™ ln 2n ~ m (l — y) for 
< m < Xn and A > 0, in which case one needs only the restriction a s ln 2 ~ A (l — y) < 1. 
That is, one must sum all the terms in ( |3.8|) below a line which makes an angle tan -1 A with 
the vertical, which implies that one must sum a large number of subleading logarithms at 
high orders in perturbation theory. 

At present, only the sum C\{y) of the first column is known, so the condition for the 
subleading QCD radiative corrections to be small is that given by eq. ( |3.11| ). To determine 
the restriction on y, we use eq. (|3.11|) in the form 

13 



as{mb) ln 2 (l - y) < 1 , (3.12) 



7T 



where we have noted that the perturbation series is really in a s /7r rather than a s . The 
condition for reliability of the QCD radiative corrections is 



l-y>e-V^. (3.13) 

For very heavy quarks, this corresponds to a region that is much larger than the smearing 
width e of the 1 / nib corrections. To see this, take the limit nib — > oo with a s at high energies 
held fixed, i.e. a s (m(,) = 67r/[(33 — 2n/) Innib/ Aqcd] (n>f = 5). Define the parameter t by 
lnm^/AQCD = t 2 - Then eq. Q3.13| ) becomes 

l- y> e -tV^ : (3.14) 

whereas e ~ e - * 2 . For large quark masses (large t), e is much smaller than the restriction 
( CT ) on 1 - y. 



As we have already noted, the residual momentum averaging procedure discussed in 
Sect. 3 can be applied to the QCD corrected free quark decay spectrum. This procedure 
yields the leading I /nib singularities to all orders in a s .f\ When radiative corrections are 
neglected, we have shown that the leading l/m^ singularities smear the decay spectrum by 
a width of order e. Thus, a reliable determination of the decay spectrum near the endpoint 
requires knowing the lowest order (in I/mi,) spectrum at least within a distance e of the 
endpoint. Eq. (|3.13| ) implies that in the nib oo limit, the lowest order spectrum with 
perturbative QCD corrections is not known in a region near the endpoint which is much 
larger than e. 



^The leading singularities arise when the derivatives in eq. ( 2.14 ) act on the Sudakov suppression 
factor R(y), not on the ^-functions as in eq. ( [2.15 ). 
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IV. NUMERICAL RESULTS AND CONCLUSIONS 



The radiative corrections become large in a region given by eq. ( p. 13 ) which is much larger 



than e in the limit — > oo. However, for a large but finite quark mass it is possible that 
3gime where e - ^ 71 "/" 8 is not much larger than e. For example, for = 4.5 GeV, 



one is m a ree 



we find a s {mb) ~ 0.2 and e - V^ r / £ * s ~ 0.02. This value of 1 — y corresponds to a smearing 
width of approximately 50 MeV, which is smaller than e. Whether this crude estimate is 
valid depends critically on the size of the coefficients of the subleading terms in the second 
and higher columns of eq. ( |3.8|) . For example, the subleading ln(l — y) and constant terms 
in G(y) (see eq. ( ^2|) ) give 



2a. 



31 \w \ 2 5 
—J Hl-y) + 7r 2 + - 



-0.1, (4.1) 



3tt 

when the electron energy is 200 MeV away from the free quark endpoint y — 1. 



Eq. ( |4.1|) indicates that it may be a good approximation to include the effects of per- 
turbative QCD corrections on the shape of the endpoint region using for dT^ ee /dy the free 
b-quark decay rate including the leading QCD double logarithms with eq. (|3.3|) . The small- 
ness of eq. ( fl.ip arises from a cancellation between the ln(l — y) and constant terms in G(y). 
Each of these separately is not particularly small. Thus we are not completely confident that 
higher order perturbative corrections are negligible. It may be possible to sum the second 
column of eq. ( |3.8|) using the methods developed in fl4| . Such a summation would provide 



useful information on the importance of higher order QCD corrections. 

The shape of the endpoint region of the electron spectrum depends on the matrix ele- 
ments (B(v) \ iD^ 1 . . . iD^ Ln \B(v)). Neubert estimates these matrix elements using a quark 
model for the B meson |7|]. Eventually, these matrix elements can be determined directly 
from experiment. For example, the same matrix elements occur in the 1/mb corrections to 



semileptonic b — > c decay and in the decay b — > 57 [15|. Thus a precise measurement of the 
electron spectrum in b — >• c semileptonic decay can be used to obtain the endpoint electron 
spectrum for b —>■ u semileptonic decay and the photon energy spectrum in b — > 57. The 
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order a s radiative corrections have also been computed for b — > 57 [[TjJ. Let x 7 = 2E 1 /rrib, 
and define 

F(y) = / -^dx, , (4.2) 

J y CLX>y 

where dT/dx^ is the inclusive photon energy spectrum in b — > sj, neglecting the strange 
quark mass. F(y) for b — > 37 and dT/dy for b — > u semileptonic decays have the same 
ln 2 (l — y) but different ln(l — y) singularities as y — > 1. However, F(y) for b — > 57 and 
dT/dy for c — » d semileptonic decays do have the same ln 2 (l — y) and ln(l — y) singularities 
as y — > 1 in the order a s radiative corrections. 

The methods in this paper and ref. for describing the endpoint region of the elec- 
tron spectrum apply when the endpoint is dominated by many states with masses of order 
^/mbAqcD- However, in the non-relativistic constituent quark model estimate of ref. IJ, the 
region beyond the B — > X c ev e endpoint is dominated by the single decay mode B — > pev e . 
If p dominance is found to hold experimentally, then the sum of the leading singularities is 
not a valid description of the endpoint in a region which is as small as the difference between 
the B — > X c ev e and B — > X u ev e endpoints. 
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